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fourth-order nite-dierence method supplemented hy special two-way arti cial boundary condition* I ABC’s! 

to sohe the NLH as a boundary value problem. Our numerical methodology allows for a direct 

comparison of the NLH aid NLS models and for an accurate quantitative assessment of the backscattered 

signal. 
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COMPUTATION OF NONLINEAR BACK3CATTER1NG USING A HIGH-ORDER 

NUMERICAL METHOD" 

d. MBIOI M. .LAN AM* V ISYNKO* * 


Almimct. I in- U'HkUiJ t i SdirurJu^r equation N’LSj 1* the i«J mode. foi | >rMpa t k.»ii of .nt ease 

laser eMins in Ken media. Du* NLS is derived fiom *li#- nonlinear Helmholtz equation (NLHi by enqnoving 
‘fa* |mj axial appioxiiua'iun and neglec'.ng the backseat tcied ua\f»s. Ill this study we Use a foiii til-order 
iiiiit fiddler ence method supplement i»\ specia* Two-way artifcc.a. boundary conditions (ADCs) to solve the 
NLU cvs a boundary value problem. Our [turner i<. a* methodob.tgy allows f«.«r a direr. comparison of the NLII 
and NLS models anti for an accurate quart*. tat .vc assessment of *.he bac-kscav.ered dgnai. 

Key word*. Kerr medium, nvc propagation. >el: foe us .1115. fourth-order method. : wo- way ADCs 

Subject classification. Apphed and Numerical Mathematics 

1. Introduction. Die piopagation of intense Ja^ei beam* (time-liciniiomc eierxoiiiagne’ic wa\esi in a 
bulk Keir medium is usually modeled In 'lie err;caJ nonlineai Srlirodiirger equa'ioii 'NLS) for the elec' in* 
tield aiup.r ude. Since l.gh* i>*ys bt*nd '«*vaid ar^*a> with libber .ndex of refraction, the non.inear dejiendence 
n f the index of rH’rac'.tii on beam intensity lias a ntlf-foc**intj effect, wherein a sufficiently .reuse User 
lieam become* narrower as it propagate*. In particular, the NLS nuclei predicts that when the input lieam 
jH/r.er < L' norm) exceeds a g.ven critic*. threshold .V . then the beam can collapse to a point at a tiuite 
propagation distance. Lor ui«m* information of >elf-focusihg. see c.g.. [O/JJ. 

As la- >aUi propagates .t induces changes in the optical properties of 'he medium. As a result, par' of 
*he incoming wave is reflected back, a phc.iomeaon referred to as buikiiut*tring. Wry j't.e is actuuJy known 
on backs ca *. ■ ering in nominear self-fc^asiug. excep: for tin/ genera* belief tha' it *s "small." Since. however, 
sinad-magm' ude mechanisms car. have a large effect in sei-focusing 0). there is a need to accurately quantify 
:he magnitude o i backsca: ering and study inkt* tins phenomenon mav affect :he beam propagation. Another 
application winch comd great]' benefit fruin better understanding of backseat .ering j> remote seus.ng *A tlie 
atmosphere 12 . v.here the measure*.! dgnai is exar.y :he backseat tered wave. 

ihe backscatrered wave js neglect ed m the NLS model winch oul\ describe* the forward-propagatnig 
wave Cax'ula'.ou of backseat ter ing requires, therefore. going back to the nonane ar Helmholtz equation 
NLH) hom which ila NLS is derived. The NLS is an c'olu .on equal ex* with die spatial coordinate in the 
direct «oii of propagation piaiiug ‘he ro.e of ~ one." Tlieref ore. lie correct mathematical formula . jju for 'he 
NLS i* 'lie Cauf ii> .nr. a. 'ulue) problem and as such, solving it nuiiier.caJ.' .* a re.a'.nely stiaiglit foi ward 
computa ionai piocedure. In coin ladi*'. net ion to that. f he NLH .■> elliptic ,n .t* nature, and a social 
mu. fid JucuM«4ral Iroiindaiy \cjue pioii.effi iieefhv t<# \x* foirnulate*! and »mi'ei] fur ibis equation, wlucli is a 
juu# h harder r^-k fi*«n rLe - aiid|M*int ««f co*uputing. The first immeric^l s.mulaf.oiis based <»n vih uig a 
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?»*.n«U«* » pr«44*m :«* > Mil *.ep* rer*ntl> je-rf nmed j t 7] iM»x an advanced fourth-order method. 

hi tha' **ud>. *h* 'ksign f<*u ?b- ird*r rate of the i:ie*b*»d was (n/rrobur^f^l exj>eiuj*i]Taliy *111 a 

model Imear problem. Subsequently. t series </ t he grid convergence 'ests wen* conducted .n the numinear 
regime. I 11 the current |*aper we go beyond grid-couvergeiite argument* and show that the asyui»totk liiait of 
‘he MH solutions obtained .n '.he sjnuiat.ous is the corresponding >LS so.ution. lliis comparison provides 
strung supjHjit ilia: the calculated MH solution is indeed the physical one. In Section 2.3 we obtaiu an 
asymptotic estimate of the magnitude of back sect ter. ng and subsequently sirow m Section 1 that it ugiees 
with he calculated vmues 

W'aeu f he .ni.a. datum is suttic.ently large tne MS dilution deve.oj** singularities at a tiinte propagation 
dist ance see Section 2.3 1 . Since, however, physical quantities do not become Jiimite. a uatura. question 
^ whether 'he correspond.ng solution of the MH exists globa.lv. Ih.s fundamental question has been 
open f many *»ears. Ihert* have ‘xmi ardicat.oas. hnigb. that so.utious to the MH exist even wlien 
tla. corresponding MS solutions h.tw srnginar. * vised on both uuurer.cul solution *A “juudJuxT MH 
« l radons 12 4| and on asyiup.otic analysis 3J. but ib«*sc studies d.d not take in o account buckscuucring 
effeo. Tlieiefon*. out long- eim goal is *0 solve iie NLH for *liose .ncofiniig signals tha* .ead to blowup 
in the MS Mod*.. In *h* rurrenr stud}. Imwever. w* cmren'rate on * lie- mure attainable goal of letter 
nulei Main ling • .n tviins *A l*>tli analvsis ami numerical simulations; the regime when the corresponding 
miIum u iA the .\LS d<i*^ 111 k li.ow up. Our liope is that this under standing wih eventually allow. to solve the 
Mil for “anv** un «<mng s.gnai. 

2. Mathematical Models. 

2.1. The ^Nonlinear Helmholtz Lquatiom A tvpica. exiierimenta* setup (both physical and 
numeric*.) for the propagation *A waves m Kerr media is shown m Figure 2 1. .Vn incoming laser beam 
" ith known characteristic* .upiugcs normally on the p.anar interface : = U between the hncar and the 
notihiicur mcd.a. Tin* electric field £ s -rp-i • * ) *n f. l ‘ is governed bv the nonlinear Helmholtz 

equation 

2.1 1 lit.. - /-f-il. <' -A, 'il -» £!*'». •/. /j. .if.?/’ 1 i>9. 

whe.** i f is the wav*num?>*i. * - ki/.t'j, it 2 .s t lie Keri coeffic^nt. and — d.. £ . ♦ . . . ♦ 0,- : , e , t . is 
-he transverse La p^ dan (the dUL'action term), see, e.g.. j 3.Sj. For >:mpl.c .ty we consider from now on the 
ey Jndricaiiy-symmet r.c Cij^e where L — Eir. r) aud r = ^/jr* + ... * rj,_ ,. 

The ii» »u lineal medjiia ••ccnpie* th* seuii*siMce ; > n u#^ Figure 21). Omsequentlv. th* NL1I <2.1) 
nas to be supplemented bv boundarv conditions at ; — (.) anti ^ -~x. We requ.re that as : — > " x. L 
lias ^ left- traveling tvm^mtmts and tl*a: tlie |iro|fagation is ihlfract .on-dominated with tla- held amplitude 
decavinz to znuo. un. max Eir . ; 1 - (J. which al^^ means Jm k 2 - kf t In otiier words, at larg^ r's 

■ x le r< x . • x 

In »Ar .';u sm^ld be a linear *u|*rposnion of r.ght-truvehng waves. Since the actual uumerica. sumnatiou 
*s carried on a truncated AniU'Ui M < : ^ l A 'Figure 2 1). rii#' desired behavior *^f the ioiu'^u as 
, x Mx s to be captured lr a far-he*d ar .hem. boundar*. Condition 'AllC ; at tj^ urtihcid b*juiMiary 
— Ihrs IxKiiddiy <1 »nd» .011 si^ild gmuaJilcc a idhAliottlcss juopagata^ti *A ail 'la - waves raveling 

ww J ; - - X. Of -Ii. boundary »*«mditiom des.guet.1 1.0 ensure lie lamparenc' <A tliei^utei k^uudar.v to 
Id' iHitgo.ng uh'.<s aiecvJ.e»l iniiuitwi* hf/nudai y ntudiUon* 

lUr s.iua'a^n is mon oxnjdrcatcd u the mterlaa* . - •) where the total la/.d L r.U) is c«^mp*>a.d 
of a g. <*r. inconung i.gn -tia'.vl.ng) courj^xr* *n( <r Mi ai^d on unknown bacist;oncred (lef.-tratchirg# 
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A* Midi. a*- 'juuijdart comniun a' : — ii La*> :o *J<ranee the idiectionies* probation of any Jeft- 
‘fa* **Im 5 v. ave tluojj^k *.be interface and at the same :jpe be able tv "urreetly prescribe the .neouuu,; i^na*. 
Iwp-t mentation of such .•!£?(' ’.vu* clone in [7*. 

I .uuJy. we assume >yiuineuy ui r =. »j and vuiuslniu' o. ihe c.mric held us r — > ~r Ju piuctice we 
unca e .lie domain «ii som** liny,** bu. Unite and ie«|jiie ilia £ : j — 0. flit* ustiiica.un for ilie 

im* 'H* h.s appr««rh to treat the -lateral" l>* * j;i«ii ri**s can !»e found in 10 and the biblio^i aphy *heie. 

Le* u* <*Jso imc that in this **ud) ap do not take info accx*nr 'lie effect of diw;ontiiiiiity in flip index 
«#f ;pfiaiti<*f. >urt»- tIjp in'erface *. — II. Thu*, we assume tliat f^ r.O) is the inem.iinig wave after it has 
<ii;ead> |jas*ed tlacmith tin* Jjterfaee ..e.. at . — 0~). We also assume “hat .eft -traveling i.e., backscaftered ) 
wave* ere not roller ted M he .nterface ; - U back into lie domain ; > 0. One ran exj»ect the latter effect 
uj le* *ual.: .t will be ixiveMigafed in a future s* id y. 


i.2. Paraxial Approxhuaiiou and the Nonlinear Skhrddiugcr Equal iuu. Let n, be the initial 
Aid'll •/ in. .tipai^n^ U ani. We fus: Jitioduec the ciuncfisioiAJcs* quantities r. and r* us 

e s — . : — — 1 — . £ - ^rjA-j-r 1 ^‘ Ir. : • . 

r #i *!#/»/* 

a here £/./ — *v, r* i> :ne d'ffruciifjn length. I lieu by dropping .he *iide, we obtain 


J.dl II - X.#* * |fj**r - -1/*'-.. . 

a he .** / — .h i,k„ *k' 1 is ih ** non/MttLnal^ly punir trier 

lie* standard delation of the .NLS ba *-*#•! <xi .la.* asMiupfu^n Unit tla* eii'e.op** #• ,s •»!•.** l.» .ar*ixi^. 
in tha* rav*/ one can iiejrjeet :lie *en;i on the rixli'-haieJ side of 2.0i [i.e.. a|»|*l> 'lie ftunui'd ap/m/rtriaf/'#/iJ 
ami *jlit«iin the ix»nliie*a. S»ii. r *j»lin^i e'|mti#ai 


J 4i ** . r|*V-d 

fla* >LS i» in supplemented the hutud ron<ii'.«*U at ; - o 

*•</• Ol - ' r -/.': 1 . 


j 



i* need* to \* urr^/tted f* a <iht< mthin. * here th»* djeitioL of pnprfjc<r..on ; 

;•!*»* a** r« 4 e »>f ••in#*. \\e ree:i.|»h**i/e that t*u km atfrnh>j rffttln air not f aA^n i n*o octinnit by thr \LS iJ.jl. 
indeed. on**#* ( 2 . 4 ) .* so.vetl. the ovpiail solution. according to » 2 . 2 ». is the *.owly varying amplitude r times 
Ih fur* art J prupaga* Ja* oscillator) com|»oiieiAt c 4 *' : . 

2.3. Preliminary Analysis of Backscatteriug. To the best of our knowledge, no accurate 
«(uant;tative analysis u< backseat tering .n noiilnrai sej-focusi ug lias ever been performed. although there is a 
geue;* u be.ief 'lur the magnitude of the backseat tered signal is small. in tin* section we present a prelun mar v 
nsvuiptuUc study of barkscaiieiing. To do that. we consider a more general ansa.z for E t lam 2.2 1 which 
is comi****] of Imtli torvaid-propugiring yud backward- propagating wuve*. i.e.. 

i.'» L - '.r.itf I"* * ^.4<r. it.'*- - 0'r. . 

wher* A ami U w 9«t*‘.v-v»rv*ng envelops. Substitution in tin* >LH «2.1< yields 


«-“• |.l 2ik, t A ±_A .1 <‘ 2 *' 0|*mJ <-*' |lf i,k„U, S.U v- :,k ’ D\ ir U 

Changing to le- leaidimensioiutl %«n.al».e- «2.2i gc.es (after dioppuig lie tildes) 

.4-. >u. + ,l, ^-u\'-D 

Lrl Us aver**** a** .ast erpia*..on over one fast (arXa’K.u. Tor example. using la' .or expansion. we obtam 

-L f ‘ll „!<*,•. -£7 / '[« :i - U - :i/l */i(:i!l *Ol/ 2 » • 

*/* i 4 '/ / -«/- 4 


-U. 


-0 . 


Simhaily. 


•1 m- * - * f* 4 

- 4 -/ <• * '" vKi</* - -77 / [.4 :J. 4 .U)-...K 

"/ 7.-,^ 1 -/* 7 * 

- —.4 <•' * ; - (J< f‘t . 

u 

O^is^iuemiy. we ov ain tne foih/wiug c*pia:.ou for lie backseat, teied w ave 

*<ii /-II. . - .71 * A./l * 1.4 - «— 1 fl Uf*U - f J r . U ( i)i 

w It**.** 
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where 



l * i.i. - ^..4 
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Siu * t h#* st A j*.na 


•n. - • .4 II. W:i\ & M 

;j» it 

«s II - *l. »he alwr.i aaaJvs> vu^v* iuj 


U'A^Uif 7 ! 

if tl*f quitch a.e */duiig L v i , i can be confirmed Uid^^udfiitiy. then .t •xl.I a Convincing 

justiiica*.on fui the <iv* jmp'.oii *iiat backseat ter ing is indeed *i.iaJ.. Hwict*. in uur sj.iulat.on* we exj#ect to 
see »ha‘ 


2.m 


«*rj*8) A - v |£ 

Ml 


A\ 


- (*/-■). 


Ila* uuueikdl resul's %A Section 4 do wrr^imrate theae expectation*. 
Tbe above analysis alv; shows *.imt 


|£|- 1 Mill - . 

I .i I 

Therefore iie auipl.tiide *jf • he NLH vjlj .on ruo ha*, e i.pples wi li the vtatejeu£th r / A'., due lo 

'Xi top of the >^« 1 ' 'i«i ' .iiy. amji.it ode of tin; foi ' , ojd-|»io|»aj^atuiK waw.*. NLH *.niyi«it*ou* 
* Jijjfrst in*, that li.s indeed nun be the r«iy* lum* l>eeii rejmrted in [ 7 ] (see a.-**.* Figure 1.2 in Sect ion 1). I* 
#* in#* c.eaj. ii«#v»e\e| . *.# v. Iia* t x>iit Ie* iij#j»Jes observed in [ 7 ] aje a iitjueihaJ artifar* due *4# plac ing the 
.‘ai-fn-ld air.ihr.al IhhiikNiv ton ilc*sp r«» the in manual seJf-for<i n^ ^ojte. Niimemvd ntudy coiidiir'ed in [ 7 j 
•lid. iu fart. Uivolve Ur analysis *»f b *r.v rh»* l»w>«f,ioii <#f in* far-fMd boundary affe** tb<* so.utiuti: 

platum this bouudar** fir.iiei tind further away catiaed the rdudiuu in the ripjde* magnitude but never 
ahowed to eliminate thrai compieteJy. This may sti.i imply, though. that tlie boundary wa* no* "sufficiently 
fa i ’ a* ay. Therefore. no dehiute eoucl\is.oii as *o the prepuce of *.iie ripples in :be NLH solir.ou* 

can be made at this tin.** and *n.> question requires a subsequent thorough study. 

2.4. Nouparaxiaiity and backseat teriug. The fradr.onal wa\ of introducing '.in* paraxial 
approximation is rejnirted ui Section 2.2 when* f ln* right-hand side of equation 2.;J) omitted and tin.* 
NLS .s d«ucd. Iia nax«* «:aieiul uiu h so o. Sr 1 ion 1 3 ^b'^ s. b»/v.i.n.*i ilia' P.w# apyrorirhatwri* .n*. in 

l.K-.ni#. made h'-n tbr NLH .> approx Jnatrd vitj* the NLS: S<-%b rUu^ .1 (the paia-ual appiavima.^u 
iu i»e uari«uv *env. # . i.e.. as .t apjilies o the foi *ajd-j>iopa^aiim, V4a^esl and hejijeei inj^ U ( backs ea erin&i. 
VV#* reraJ. 'ha' jireviou* st idi*-s [l.#J Mj^»-sre»i that n«#i.pai axial in of tie* r.Jit-'iau.ijiK <i.e.. . 4 .. in 

* 1 »* %ense »»f S*« r wm. 2 ,. 3 i ai r^sfs *iie i olbi|#se of fin* NLS solutions. i#nt tIk^t sfudies did not rak»* into ar(r<mi ? 
r#rt< ks4vitteiinK »df#*rtH. Having said tl^i* . we stiil u*^#* tjia*. the *epaia:.t:u into nonj)^tr.<xnil mid barkMvttterjni 
tiff's, vn.rh is bay*»: on the an»utz *2.u». is vumev.ha* artiht^l. suee the probleu; in n'unluieox. Therefore, 
uhen M*e compare the NLH and NLS s«alu*.ons. It i> not preaseiy clear wiiidi j#ar» *A the d.ffeience comes from 
noi.par axial effect * for the n^ht-traveiin^ waves, and winch one from backseat >111115. A iiutablt wept ion if. 
fuju<:n r. at : — U. Mu r*. *fr< dtjji nnct ortu*o:h SLH and SLS whiti'w* i> *0(1 1 y due to tochfcatti'riuy. 

2.5. Critical NLS. i* is veil kntuvn rju-.r yumtbais of rln* NLS ' 2.4 can bK*»ui:ie s.n^u.ai when either 

/; 1)^2 tb # * Mjr.'ftntica! NLS *ur when <r U ii — 2 *ne ttiUml NLS <I> .s tlie sj>a#:e dimension'. 

H«jv'«*'oi. wn#.T‘*as in suj##.*rcii'icai colhipsi noiJ.n'uui * <hmnna .«■> o*ei diibac.ion U‘*ui the ‘'.n^umii'*. in 
Jn rir^-al coha^r-* noiJ.mvn > and diJiac a.u are al;n^ luaUuitrd win ih" sin&uLiri.y. f.V'ns'-qu* n'J' Ua.* 



smjjijiM Uun'i »;i js b.yjii- srstf.tivr to small |**r jr*Mti in* in the tr.Uml rase. but mud less mi ^ the 
• n » a*e. 

Ib** phv%*r,d ri»*f forre*j»on<W ••» the pn^a^'d'** if .user bfvua* in bnlk Kerr media i* *lje cri'.rai 
•w. as l) _ .$ and ^ — 1. However. in older to ihIu(t t|*e nnj*plex.f ’ of the 'iiai|iutH*j^. Mmv ue consider 
*bf critical ta*e D - 2 and (7 - 2. Thus. the Mil for £ - Hr. ; i and the NLS for r - r r. ;). which un- 
solved numerically ui this stud; are 

2‘ h £: -£Vr^frjli-#|£| 4 )£*U. 

and 

2-IUi /» • . -*■ r r . *e |» -| 4 # • — () . 

res|>ecuveiy. 

3 . Numerical Methods. Tire NLIi 2/Ji .s so! red us .14c ftwTL-urder Unit* difference*. Tlie choice 
*»f r lubr^ir'l'T method is piunarily by tin ne«»*?iMfv to resolve a small-scale plrcmifiieiioii 

backseat ter mx; at tL<* background of tin* forward propagating waves. Tbe NLS (2.1U) ^ also solved by 
c founlr-order scheme: it ls datura* to expect tba*. this v. ill leave .ess* loom for potential purely numerical 
discrepairr.es between tbe two techniques and as such, wii allow for a more accurate comparison. UeMdes. 
if is general.? known that a.gh*T-order methods |>r«>\*icl*- for a better resolu'.ou of waves. 

3 . 1 . Numerical Integration of tbe NLH. Our numeric*. morbid for solving the NLII j» delineated in 
[ 7 j: here we u*dy outline its key eminent s. We use a convent ional fourth-order cent r aJ-differeuce di&cretization 
of tbe Ldpitff.iau: .n so domg tb»* stencil is five-node wide iu both r and ; directions. .Vs the equation is 
nonlinear. *.h impleiiteu* a nested iteration scheme. On tlie outer loop, we freeze tiie nou.iaear.tv. i 
consider ire coetfir.eut k* *A <2.1 j a* a given fuar.wou of r and wlacir is actuary obtained by faking 
£ 1 fi din the prev.ou* i'erat.ou. Ain* wav we arrive at a .inear equat.ou w.th variab.e coefficients. Tlie 
.a* er is also suhed bv .feint ions oa be jjier of tl^e nested scheme. Heie, we leave tbe eutire \*arvai^ 
j»*j > ol ia 'sp.a .a.»n. wLkIi js pidj»«.<Uotui. «o <. on tin.* .own x*ni aind on la: upper .evel need o iuvei . 
onlv i U cons aja-cvlbea:n». .aieai Hclmla^lu (operator k^J . Formal h. «.*ui iteiauoa x/beua. reseiiibks 
la* bxHd*jHi.ui i*ppinacii !iov e’ ei . no rigorous con* eii^euce iliet^i\ i*» a^a.lable \et. awl the cunt ergence 
nas to be a-*«.*w<| exj»»fin.mtaih . *Di»* advantage of isii^ *iiev* n»-s'e<J .relations is that firs' . t|ip i i#.*' Ii^J 
e\nr i/dl'. iHilia#-. ’ti »lie j»ia*^.t**d M»|tr;on «if one .md tla* h;iiii#» liia*#n constant coeffiiicir e<|tiatioii c!ri\»*n 
b;. dJfereu* ** >iuv As expbiiiasl bel*iv. . tjjs c an be done effic.mth on *be dixiefe U*ve # . Setund, the 
radiation Imundar;. and tla* tw<>-ua> AD^'s aie moM convenient to set cm the upjiei r.ua* level <»f 

b* .feint*. m xbea.e a.re/id*. f<n fb # * bn ## ai »-4iastajit-coefhnent «fpeivitor. 

To s* Jve f be Ai,ear-c his* m.f c«cefbc.ej|t lleb.jlaii'4 **cj Mfi m 'dix refe / (Miifeipar* of - //. where 

y is the i ii^jr >hdiid side generated b<* previous .terationi we fir st separate tJa* var.ab.es by implementing 
la dixTete Four.er trai isbxu in t|*#* transverse clirectiou /*: tb" boundar coialitams are syiainetr> at r — b 
diiii z»T'.* Dir.diier at r — r, J< A . in.s yields a rolas '.i.m of fouitb-<Xfler oiie-d.iaeus.oua. Huite-dJlerence 
ec| jafi*xis parame erued b« la. duai Fourier tar-ab^*: each of »la* .att^T needs to be sohed independeiitb. 
lbe wo*wav and radiatioa AiiCs a' : - U and ; - respec ^eiv. are x*' in the Four.er sjiatt* as v.el., 
ijr. -»'j»aia W*. foi eacn *>»( * b e ab>renreirti<ji*ed ^if*-dimensioiuil e'juaf ions. Ibis i> done lr. hrx ideu'Jvini^ 
b< l.wvjh -ueltp uaivip <*ijy:a'iU'>«b> foi »he b^aoip *ia. ous vuxon of each <ji4f>dmieusioiiai (.sjnu.a^n. 1 
.> m 4 eii «i.f o not I i/U r • *ii ia^u^b 'la «»ri^nuJ diiluniUaJ fjLia’am .s oi be y.t'ond ordci. v*e jrr 



i % .ug **• app«‘Quma* # on and as *j«h. each hfJMograeoj* di j*;etc one- dimensional ef f ua**ou 

iir*A f«Hlt llMMfh iJldeje'ialenT NMJfliUiV Oil#* pa*T A 1 Hlb laft ej i* ppl o.\j!j,a’<** ' h* genuine IlilfflH <#f the 
drlf'T»n'.a. ^pitr. on. tlm uiay be erher uavelmg or evanescent wave* depending on the value of the 
dual Fourier variable. 1m. o*hei pair it a pure* nuiaerica* artifact. these warn arc a] wav* evanescent. bu; 
‘heir presence implies t|*at every t* tt* equation ietiuire> *.wu more boundary conditions compared to the 
original differential equation. Iiu iad.tr ion boundary conditions are constructed by requ.riug that on tine 
*eft l»»undarv — Uonly the left -traveling and /or left-decaying (evanescent ! waves be present in the solution, 
and on he r.gh* buuudarv r - only the right-: raveling and/or right -evanescent wave* be present in the 

>ohj .on lid* election .s rendered bv » la. vxca.led one-way discrete Helmholtz i.t |Uitt ions which ore 1 la* lincur 
iaxaogt'neous i**Jaiiou» ‘hat 'Min** 'b»* span of all appropiiaie mode* for era.h boundary. The iuo*a> ADC 
hat a.su prcsn.lfes tlie incoming sj^nal at . — d b constructed on the luasis <jf the euiresjiondiiig radiation 
ooundais condition In Mibstr'iting tin* light** raveling incoming wave into tlie une-way-*«>thedef Helmholtz 
*«]uafi<ifi and as Midi creating 'Ik* .rihom(»gen«*rv i»f a particular form, m*- [7J. Simple consider at :«ai* based 
oii tin* Lnear super pmi* am pr.ndp.e and iimqiieness guarantee diat the resu.ring iioninJinogenintis relation 
a ill correr*./ speedy the urconnng signal a: ; - (l and still ensure tire reflect ionics proiJogat.on «if all the 
outgoing ii.e.. left-traveLngl wave* *lnough : - 0 *oward ; - - x. 

Ao concerns *he computational complex, ty of the resulting a.goiithm. if we introduce the grid dimensions 
.V and .V.. then the cos*, of both the direct and inverse FI T wU be GMA'.AV iu.V r b The coat of solving 
each *A the A, one- dimensional systems wiJ Ik* Lneur w.th res]**** to A.. Indeed. in the course ol iterations 
each of these svs*ems needs to solved man;. *.mes for d.tferen* right-hand sides. Conscqucn*.y *he sjwrse 
LI. de<.*om|»osltior can be performed only once ai^ad of tune, and the cost of each backward substitution is 
linear AJ'ogrt her. the cumplcxry of each it era* .011 is s* J. C/t A A r hi Ar). 

3.2. Numerical Integration of the NLS. The ALS i^.lDt b diMieti/ed n tlie r direction using 
stajuUrd fo irtli-oider centra, differetjee*. It is iiregj at«*d in . w.»ha four-stage Him ge-Kutt a me* hod staring 
with the inJia] data L Uk (r.Oi. The IxHindary condition at the rename la ter a. boundary r — r tUl , K is zero 
JL).ricb.e* . as in the case A the ALII. 

4. Computational Heauits. Li. ace»xdaiK*e with he d.scu»;on of Sect^KJ i we have designed a se* of 
nuni^ica. experiments ajmed a' a«inev.ng two object n*es: <i» V alidate the computational aJg«Krithm *A [7 # for 
solving tlie ALH iiT'^ugh a cohiimrisuu v.ith numerical sjnulatioiis A tne ALS model, and (II: corroborate 
hc< 'xicksciiuci Jigelli cts cap ixed la the ALH UKxii l sc*a.iv|'Jiairata:ai.v wi.h the mx.pui axial. ty puiamcicr 
/ as suggested In the auai'sis o. Section** 2.3 l 4 Kegaid.ng the hist tlijw.iivt'. le us note hat previously 
ae ha-.e eMtsI he nnmeiicai aJgoi.rhi.i of [7] in the nunhnestr regime using grid ctimergr^H e. but never 
cunipajed .t v.ifh au\ «ab#*i ajgor.rhm for computing *he pjojiaga'ioii of waves in Kerr media. As concerns 
he MHiifid *i4»je«tive. .t amounts to 'lie accurate numerical computation of backsca* rering in nominear self- 
fiH iising. njid w# aje t uriently unawaie of any pi evh hi s t«hni(| ie with •'U.iXir capabJities. 

To be j'h** to coaduc*. an accuia'.e coniparxm tiie numerical predict.ons tamed with the ALII and 
AL*5 :uo*J« ! te*iuation> 2 f )i and «2.1M). respectively i. we have chosen a regime with the input piuw#*i b*:.^' 
critical, for which the so.uuoi. of the ALS d<>es m^t devc.op singularities. Ue take l*o - * oral * — 0 U4. which 
corresponds *o 74'/< of fin* critical je.»wi i A . sec 7. 11J. The incoming beam proli.e is L u . t <r. Oi - c“ ,r r "‘ . 
* jc:i tl*a* the l>e<ii;i width r (J iv mjt.'h le»» *han e lll( , A . lo alloy* f./r th“ vanatam of tin* nonpar uxiahtv 
paranie er / — *dvy. we 'ary *ne »jeam width r , whhe jn***j»iiig 'lie wa'euumber Ar«j arm T he f|uaJitj v eoyT - 
h.J ojl io.s le* liuf .oua. critical jh>wi i jiRhaiigfsi. 

In FigJJf 4 1 ,,vf ‘ sln/v ht on-a.'-.s ampi.'ud'* pioILcs i ^i h** ALII and ALS numeric a. su.utioirs. 
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Figure 4 1 leiti corresponds o / - i/s and Figure ill right) corrojoud* to / — 1/10. We ji.ot the 
vu.u<*s *4 tin* computed solu r -ou '.ha the uxa> symmetry r — U because tin* is the mov. iirc*rc»tnig locuion 
;u ih* 4 doma.ii nheie the geuulnel* uouluieai phenomena take place. A dear manifestation of nctu.mear 
•*lf-fuc j*.iig i* »lie ~buiup. or |*eah on iie sohit.on nines ill Figure I I vhose »aiue is hHjhn than tha* uf 
*Jh- .nn * 1.11114 wa’.e .'O.'Ji. 

I* s *■» w the! • • *r forth f — 1/H arid / — 1/W» *hr VI M #*nf! VI.A iuive« ; n Figure l.l are d<i*e to 
«*ne a;e»r|n*i . As the ALS is a iveil-esVtlih sited tin nip) r has regular so.iifxms for MibcrrhvtJ ill.tiaJ jwneers, 
we mne.ude tiiat out numeric al aig*Hri iiiii fin solving the ALII 17) [that starts the .teratiou process wi'h tlie 
.in f id. gjes-> L — U indeed converges to -.he toirer. k4ut.ou 

WV diso no*.ce that the discrepancy bet ween '.he ALH an<J ALS curves on Figure 4.1 left) is larger 
hai. tha* «.*n Figure 4.1 r.gh T i. Iiu> beiuiv.cn .s exj*ected according t^ the anaJys.s ol Sect.ou 2.J. because 
he disci epanr I.h am u the T, *.o cur^*> i* due to nonpar axial* tv and backseat ter mg. in particular. these 
snnuiarioxis sjgg^t tha’ ae .VLS is indeed the av Urjrfotic .unit </f tiie ALH as / — - ♦ U. IVibajo the tiros: 
apparent manifestation ..f 'he pr^nee «.4 ImckH'tif terjig .n the •< th>* ,\L11 is that h*- ojcaputed 

\a.ue of 1 li » 4 i «> ul ejec.ru* held d< . — U dijjivi flow ha« *^f he iu«*oiiiiiig wave « 9 one can dearly m 
F* ^ur* 4 1! ahr-re »*< nxmi ua <>u tjn rvo mrvis ub'aixicti lor / — 1/S TIac* siaai. np|»ln .n 'lie ALH s*^l j*aoii 

aia.v a.v.» fje e« i deuce A bad.»“*a' mag (six* S»* <, 'i«ra 2.3|. 

Ae A t. de f|uaa ifv t h#* oach*n «it ter.i.g eKer* h*. rouij»iif nig a sene><if Njiiitkdi jiaii-* i ALS and ALHi for 
rtih:it.o;i/d va.uenof / In Figuie l..'hl»d' I we sl*»»v Ir. as’etiskson in* log-kig m-jm* the <{iAautj\ 

4 1. [<k~r j y ’£ U Ml niU Ui _ 'll 

cf foiiau^i <2 S;;. where £* »he cxnj*|»tit e*i t-A&juu ol tue ALH and e .s tjie er.ut knr to the ALS tlia- 
tae jiitia. eond;tioi. imO.U) — 1 Uh / — 1/10. 1/12. 1/10. 1/S ihe *ohd litre on iigrre l.^lef* ) 

h«.t ti* . e.i.»y*h tL<« '*iriirjeited data .s •! ~j f iv \ lias e>s^|/ uili' - cun»^H^n:^ f Iatt *.he inagiii* rde of 
o^-k^vxtf e.u.»5 iird**e.J v/di'-s •jua*]: a tica.lv v.rn f Jn* nonpar axial.". jni/aui#*“i /. a> preuicteo in Sectwn. 2 d 
u.<i .ur'ra* d*v*ussMj nr Vrtrow 2.1 

L* js alvr 'Ira’ o *.a-uae ia. jru.u» < I . ) ue do ir*.* rea.l. ii'«d to solve On ALS bnausi* t In* 
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iiu’ia. ** ; - 0 1“ jgv»*a To tLi* actual r*»Ui|»iit*«l solution* of tL** >L11 and >LS. u#* j#M on 

ITguii* 4 . 3 (ri*ht i tb«* (| joi^tity 

47i ‘|£ m :) -|i , M.ril| 

-«.*t / — i/> ami / — i /Mi (tttnjt* *<»1 jp* as «.*t4 F.y.uiv 4.1 1 as a Iuju .on of .Jip normal. zasJ I’l opa^a *ou 
•iistajirp f« »i . > n. AHhh^rIi tin* rum* on Fi^wt* |.3(ri&lifi an* uKi.lau«ry. up still s«*i* tlia’ that tin* 
•Jifhicnrp l*?* wen *1 k* *olutiut^ of Hip NLH and NI*S d«u**a%ps foi all ; v.rh flip d»*T«»v7 of /. 
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0. CoucIumou*. U* Lin* o/uipuiod ujxiipj j‘ # nil\ Ik >o.uiioias to 1 li« * imaluvui S.iii • ikiui£<. i •‘•juu aoii 
non* in* aj H ; .uJi* *li z ptpjauon ho n *jf a*o<l'*l la. jiiOjiiiif,uUon of 'in* 1 ‘•Lai mom*.' « lro i un*ii,'ni/< i* 





**''•*• »i. K 1 *!! I Mil a .is ^ »/ hi u*.i»k (i u** r - tmrth. u h me*h«*i Mipp.e minted If he tai*»*#i' 

temiKimv ' • indr.Mi* tli*' KWiranfe** 'll'* |»r«v**r «>f f|*e <•» fl^y entei <unl I*’«\t the 

uttuputaMomU domain. A» ie* MS A ;ott»/ieral an e*tabh>h*d naMel. *1*1- a*i cement the MM and MS 
'Xiuldtiou* pr« r ade> a juMifiur.oii that the MH a«&orithui indeed converge* to the correct |>h'> 4 cuJ 
Hdiit.ui. On lie other hand, he MM is a more compreheitoive mode. tlu»t. uuhke the MS. take* .ntu 
en .'iai tie* phenomenon of lejfilineaj- harkM*a**erin£. A* >urh. ne atiJmre the Mnah di*nepaj|rie» that do 
U*r »ve*ii tn»* MJI cual MS solj‘.un to nonpar axial ami bac4*oiuer.ng etfect* U\ anaJ'^iiig several 
eompu* a’looai vai.an‘9 that correspond to ddferen* vtnue* *A the uofiparu.\i<d.^ iwrame’er / ve have beeu 
uhl* vouo e.n i • e Imt che um^u'.cidc *.*f the iMi'.iy.liUtit.d v av* 1 .ndtvd vale** ‘pmdiu .rulA wi li i Li*« 
j#aj aim ei *eeuidm>., o tlm tbeot't ical piedicamv Io the beM t.4 oui hno*h*lgc. thi* i* •In* luv stud* e%er 
hat n.k/* •> foi an HfT'Jictu* fpiiur.ita* Ae «Mima:iun «jf l«aeh%a "etm^ .n uouhueai te.f-focusing. 

Ackno* IferigiciPtit. (». Filmii u-riJd .ike ti» 41 km** l#*dm useful iii*ruvh<m* t*. if li (». Papain* Oiitotj. 
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